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Vortex lattices in three-component Bose-Einstein condensates under rotation: 
simulating colorful vortex lattices in a color superconductor 
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We study vortex lattices in three-component BECs under rotation, where three kinds of fractional 
vortices winding one of three components are present. Unlike the cases of two-component BECs 
where the phases of square and triangular lattices are present depending on the intercomponent 
coupling constant and the rotation speed, we find triangular ordered "colorful" vortex lattices where 
three kind of fractional vortices are placed in order without defects, in all parameter region where 
the inter-component coupling g' is less than the intra-component coupling g. When g > g' on the 
other hand, we find the phase separation; In a region where one component is present, the other 
two components must vanish, where we find ghost vortices in these two components whose positions 
are separated. In the boundary g = g' , the accidental [/(3) symmetry is present, in which case two 
vortices in different components are close to each other in some regions. 
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I. INTRODUCTION 

One of growing topics in condensed matter physics, 
high energy physics and astrophysics are multi- 
component condensations. In this sense, interesting fields 
of research are exotic superconductors, multi-component 
or spinor Bose-Einstein condensates (BECs) of ultra-cold 
atomic gases, superfluid '^He, exciton-polariton conden- 
sates, nonlinear optics, hadronic matter composed of neu- 
tron and proton Cooper pairs and quark matter com- 
posed of di-quark condensates consisting of quark Cooper 
pairs. In these systems, fractional vortices can be created 
by rotating superfluids or BECs or by applying magnetic 
field on superconductors. Fractional vortices are char- 
acterized by rational or fractional quantized circulations 
for superfluids or BECs and fluxes for superconductors, 
as found in various systems: superfluid ^He |TJ[5], p-wave 
superconductors [TJ |3HS] , multi-gap superconductors |51- 
Hj, spinor BECs ^101 HI], multi-component BECs piHSi] . 
exciton-polariton condensates |25} |26] . nonlinear optics 
|27j . and color superconductors as quark matter J28fl34| . 

Theoretical and experimental investigations of the 
properties of vortices are easily done in BECs of ultra- 
cold atomic gases [35 , which reveal to be the ideal phys- 
ical setting among the others for this purpose. In the 
theory side, BECs can be quantitatively well described 
in the mean-fleld theory by using the Gross-Pitaevskii 
(GP) equation. Experimentally, BECs are quite flexi- 
ble and controllable systems being the atomic interaction 
tunable through a Feshbach resonance [36, Also, the con- 
densates can be visualized directly by means of optical 
techniques. Thus far, two-component BECs have been 
realized by using the mixture of atoms with two hyperfine 
states of ^^Rb [37] or the mixture of two different species 
of atoms |38fHn] . In Refs. P^1 - [T^ . the phase diagram of 
the vortex lattice forming in two-component BECs was 
studied, and a rich variety of lattices was found. When 
the inter-component coupling g' is increased, different 



configurations from an Abrikosov's triangular lattice of 
fractional vortices to a square lattice of fractional vortices 
have been obtained [T31 - fr5] when g' < g. 

On the other hand, particularly interesting cases in 
high energy physics and astrophysics include quark mat- 
ter composed of di-quark condensates consisting of quark 
Cooper pairs [41J, which may be realized in the core of 
neutron stars. When up, down, strange quarks partic- 
ipate in forming Cooper pairs, the order parameter A 
becomes a 3 x 3 matrix, on which color S'[/(3) and fiavor 
SU{3) symmetries act from left and right, respectively. 
The Ginzburg-Landau description is available in the per- 
turbation theory of quantum chromo dynamics. In the 
ground state, for which the condensate can be considered 
to be in a diagonal form in color and flavor without loss 
of generality, the color and flavor symmetries are spon- 
taneously broken down into a diagonal combination of 
them. For this reason, this phase is called the color-flavor 
locked phase. Since color symmetry is spontaneously 
broken and completely screened, quark matter in this 
phase is a color superconductor. At the same time, it is 
a superfluid since the baryon U{1) symmetry is sponta- 
neously broken. Consequently, vortices appearing in this 
phase carry both the quantized circulations and quan- 
tized color magnetic fluxes [28-31,. They are referred 
to as non-Abelian since the color fluxes they carry are 
non-Abelian magnetic fields. The circulation which they 
carry is quantized to be 1/3 of the unit circulation, due 
to three components having non-zero value in the ground 
state. The interaction between non-Abelian vortices was 
calculated in (5^ to be 1/3 of that between integer quan- 
tized superfiuid vortices, which implies the formation of 
a colorful vortex lattice when the color superconductor is 
rotating as in conventional superfluids. The phenomenon 
for which the vortex lattice behaves as a polarizer when 
light passes through it [32 is one of the interesting conse- 
quences of the aforementioned properties of quark matter 
in the color-flavor locked phase. Another interesting con- 



sequence is the existence of a colorful boojum [M] at the 
interface between a color superconductor and the confin- 
ing phase in which all color degrees of freedom must be 
canceled out. In the latter, three color magnetic fluxes, 
say red, blue and green, with 1/3 circulations are com- 
bined together to a purely superfluid colorless vortex with 
the unit circulation, which can penetrate into the confin- 
ing phase. Even being the colorful vortex lattice one of 
the most important signal whether the color-fiavor locked 
phase is realized or not in the core of neutron stars rotat- 
ing very rapidly, there has been no quantitative studies of 
colorful vortex lattices thus far, because of large degrees 
of freedom of non-Abelian SU{3) gauge fields. 



The purpose of this paper is to propose simulating cer- 
tain aspects of colorful vortex lattices in a rotating color 
superconductor by using a BEC of ultracold atomic gases. 
We create 1/3 quantized vortices with three different "col- 
ors" by preparing a rotating three component BEC. Be- 
cause of our purpose to simulate a color superconductor, 
we consider the symmetric case where all three intra- 
component (inter-component) couplings are the same 
511 = 922 = 533 = 5 (512 = 523 = 531 = 5') and all chem- 
ical potentials and masses are the same as the case of the 
mixture of atoms with different hyperfine states. Com- 
bining a set of three different vortices together results 
in an integer vortex without a "color". These features 
are the same with 1/3 quantized vortices in the color su- 
perconductor. We solve the GP equation for a rotating 
three component BEC subject to a harmonic trapping 
potential by an imaginary time propagation, and find 
"ordered" Abrikosov's triangular lattices in all parame- 
ter ranges oi g' < g. Here, the term "ordered" implies 
that each of the three different fractional vortices wind- 
ing around the three different components constitutes an 
Abrikosov's triangular lattice without defects, and that 
the total configuration is also an Abrikosov's triangu- 
lar lattice. In other words, we did not find any defects 
or displacements of vortices with different colors. This 
situation is in contrast to the cases of two-component 
BECs |13tll9] in which there exist an Abrikosov's trian- 
gular lattice and a square lattice of fractional vortices in 
g' < g, depending on the inter-component coupling g' as 
denoted above. When g' — g the zeros of the density 
of the different components can be overlapped. How- 
ever, the comparison between plots of the phase of the 
order parameter of each component shows that the vor- 
tices are not coincident. In the parameter range g < g' 
of the phase separation, we find that vortex sheets are 
present like the case of two component BECs [^. In 
a region where one component is present, the other two 
components are absent because of the phase separation. 
In that region, ghost vortices of these two components 
are present and they are separated, even if there is no 
symmetric ordering as found for g' < g. 



II. GROSS-PITAEVSKII ENERGY 
FUNCTIONAL 

The energy functional of the Gross-Pitaevskii equa- 
tions for the rotating BECs subject to a trapping har- 
monic potential can be written as: 
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where the derivatives include (r, 9) coordinates. We mea- 
sure distances and energies in terms of & = ^Jhjrau) and 
hio respectively, where m is the mass of the atoms and w 
is the frequency of the trapping harmonic potential. As 
denoted in introduction, we consider the symmetric case 



511 = 522 = .933 = 5, 
TOi = ni2 = m^ = m 



512 = .923 = 531 = 5', 

Ml = M2 = At3 = M- (2) 



When g 7^ g' , the system is ^3 symmetric under ex- 
changes '^ i -^ ^j. When g = g', the system is C/(3) 
symmetric ^^ -J> % = W'i'j with U € C/(3). 



III. VORTEX LATTICES 

We minimize the free energy (ITl by the non-linear con- 
jugate gradient method (the imaginary time propaga- 
tion) in the FreeFemH — \- package. In all the numerical 
simulations we take g — 1, il = 0.96 and /i = 20. We 
calculate the ground state of the system by minimizing 
the energy functional (fTl. 

We first start as the initial configuration from the 
ground state of a non-rotating BEC where no vortices 
are present. We fix the rotation speed and let the vor- 
tices form in response to the rotation. We find that the 
vortices are nucleated as integer vortices disposed in an 
Abrikosov lattice for each value of g' . This implies that 
an integer vortex composed of three different fractional 
vortices, \E'i ~ e*^ with the other two having no winding, 
is metastable. We show some of the results we obtained 
in Fig. IT] The defect in the lattice of Fig. fll (c) appears 
by chance and it is metastable. 

Next, we introduce a perturbation in the lattice of in- 
teger vortices by splitting the one near the center of the 
cloud into a set of three different fractional vortices, with 
very small spacing. We use this perturbed configuration 
as the initial configuration in the next step and then ob- 
tain the final configuration for each value of g' . We can 
expect that in a real situation the fiuctuations of the 
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FIG. 1: (color online) Some of the results obtained when the initial configuration used in the calculation is the same as 
the ground state of a non-rotating BEG. The top pictures in each panel are plots of the density profile for each component, 
n-i = l^ip. The bottom left picture is the full density profile n = |*iP + |*2p + |*3p, while the bottom right picture is 
the same plot, where the various components are distinguished by color (gray tones). Red regions correspond to zeros in the 
density of the first component, blue regions stay where the density of the second component vanish, while green areas are the 
zeros of the third component density. However, no different colors (gray tones) are visible in this picture because vortices of 
all three components are coincident. The value of S is reported on top of each panel. In all cases we obtain an integer vortex 
lattice. 



trapped condensate will be responsible for the pertur- 
bation needed to eventually break the nietastable inte- 
ger vortices into the fractional ones. This behavior and 
the final structure of the three component lattice depend 
on the sign and on the magnitude of g', the coupling 
constant of the inter-component interactions. To char- 



acterize the system, let us define the ratio of the inter- 
component and intra-component couplings by 



S = g'/g- 



(3) 



We examined a wide range of values for 6, from Sn^in = 
-0.3 to (5max = 1-5; with the separation A6 = 0.1. In 
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FIG. 2: (color online) Integer vortex lattice obtained when 5 < 0. The top pictures in each panel are plots of the single 
component density nt = I'&il^. The bottom left picture is the full density profile n = l^&ip + |^2p + |^3p, while the bottom 
right picture is the same plot, where the various components are distinguished by color (gray tones). Red regions correspond 
to zeros in the density of the first component, blue regions stay where the density of the second component vanish, while green 
areas are the zeros of the third component density. However, no different colors (gray tones) are visible in this picture because 
vortices of all three components are coincident. Because of the attracting interaction, integer vortices do not split into fractional 
vortices. 



Fig. [2] we show the results we obtained for S < 0. We 
obtained an Abrikosov lattice of integer vortices, as can 
be seen form the bottom pictures in each panel. When 
i5 < 0, in fact, the interaction between the vortices be- 
longing to different components is attractive [22^ and, as 
expected, only integer vortices are formed and they ar- 
range in an Abrikosov triangular lattice. 

When (5 = each component is decoupled from the 
others and vortices of different components organize in 
an Abrikosov lattice independently. 

The repulsive inter-component interaction present for 
< (5 < 1 [22 makes energetically convenient for the sys- 
tem to form fractional vortices instead of integer ones. 
The results obtained for this range of values for S are 
reported in Fig. l3ja)-(c). The schematic picture of the 
whole vortex lattice structure is drawn in Fig. IsFd). As 
can be seen from these figures, the fractional vortex lat- 
tice does not show any defect and both the total and 
individual lattices have the Abrikosov structure. This 
fact states the robustness of the shape of the fractional 
vortex lattice in a three component BEC. One can find 
that vortices in each component constitute the Abrikosov 
lattice. Moreover, the repulsive force among vortices in 
different components determines relative positions of the 
Abrikosov lattice of each component, and consequently 
the total configuration is also in the form of the the 
Abrikosov lattice. This is possible because the number 
of components is three. In fact, in two component BECs, 



there are two kinds of the vortex lattice structures de- 
pending on the inter-component coupling, that is, the 
Abrikosov's triangular lattice and square lattices. On 
the other hand, in our case of three component BEC, the 
Abrikosov lattice structure is robust in all ranges of the 
inter-component couplings for < 5 < 1. 

If (5 = 1, the symmetry of the system is enhanced to 
a C/(3) symmetry. A result is shown in Fig. Qa). The 
situation is between the ordered Abrikosov lattice for S < 
1 and the vortex sheets for (5 > 1. A new feature can be 
seen in this case: density profiles of different components 
can have overlapped (almost) zeros. The right plot in 
the second row in Fig. BFa) shows a superposition of the 
density profiles of each component, distinguished by red, 
blue and green colors. However, other colors are present, 
which are generated by the mixing of the original ones: 
yellow comes from the mixing of red and green, purple is 
the result of the mixing of red and blue, while cyan arises 
from the mixing of blue and yellow. In order to specify 
the positions of vortices, we can look at the plot of the 
phase of each of the order parameters, defined as (pi = 
—ilog(^!i/\'^i\). By using a schematic representation, 
reported in the last row of Fig. EFa) , we can compare the 
positions of the vortices in each component and can see 
that vortices are not coincident, even if the (almost) zeros 
of the density are overlapped. This happens because two 
vortices with different components are close to each other 
in some place, although they never coincide. We also 
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FIG. 3: (color online) (a)-(c) 1/3 quantized vortex lattice for g < g' . The top pictures in each panel are plots of the single 
component density Ui = I'l'ip. The bottom left picture is the full density profile n = |*l/ip + |^2p + |^3p, while the bottom 
right picture is the same plot, where the various components are distinguished by color (gray tones). Red regions correspond 
to zeros in the density of the first component, blue regions stay where the density of the second component vanish, while green 
areas are the zeros of the third component density. By this plot it can be seen that vortices in different components are well 
separated for this range of parameters, (d) A schematic picture of the 1/3 quantized vortex lattice for g < g' , where the 
different shapes and colors distinguish the different components. 



notice some symmetric structure in the lattice. 

When S > 1, the phase separation occurs and vortex 
sheets appear as in Fig. W[h). This resembles the cases 
of two-component BECs, but new features arise. In this 
case, in a region where one component is present, the 
other two components are absent because of the phase 
separation so that zeros of the density profiles of the 



other two components are completely overlapped, as can 
be seen from the complete mixing of colors in the right- 
most picture of the second row of Fig. E[b). By using the 
phase plot as done for the 6 — 1 case, we observe that 
in a region where one component is nonzero and the am- 
plitudes of the other two components are zero, there are 
ghost vortices of these two components whose positions 
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FIG. 4: (color online) (a) 1/3 quantized vortex lattice for g = g' with the (7(3) symmetry, (b) vortex sheets for g > g' . 
From the top in each column we show: plots of the single component density Ui = |^ip; plots of the full density profile 
n = l^ip + i^2p + l^al^, where in the rightmost picture the various components are distinguished by color (gray tones), 
with red corresponding to zeros in the density of the first component, blue staying where the density of the second component 
vanish, green corresponding to the zeros of the third component density; a plot of the phase of the order parameter of each 
component, (j)i = — ilog(*I'i/|^i|) (black represents zero and white does 27r); a schematic drawing of the positions of vortices 
distinguished by colors and shapes, with the dashed lines in the right column being the domain walls separating sheets of 
different components. Because of the superposition of zeros of the density explained in the text, other colors appear, due to 
the mixing of the original colors: yellow, cyan and purple correspond to red-green, blue-green and red-blue mixing respectively. 



are not coincident. 



IV. SUMMARY AND DISCUSSION 



In conclusion we have found a robustness of the or- 
dered Abrikosov's triangular lattices of fractional vortices 



in three-component BECs in a large parameter region, 
where the inter-component coupHng is less than the intra- 
component coupling g' < g. Our finding is in contrast 
to the cases of two-component BECs, where triangular 
and square lattice are present depending on the inter- 
component coupling constant and the rotation speed. In 
the range g < g'{5 < 1), we find that (almost) zeros of the 
density of different components can coincide, even if the 
vortices lying inside these regions are separated. When 
g = g'{d = 1), two vortices in different components are 
close in some places, where we see the overlap of almost 
zeros of these two components. When g < g'{6 > 1), vor- 
tex sheets appear similarly to the case of two component 
BECs. In a region where one component is present, the 
other two components must vanish because of the phase 
separation, and in that region there are ghost vortices in 
these two components whose positions do not coincide. 

Our results imply that in the region < (5 < 1 trian- 
gular structures of colorful vortex lattices in three com- 
ponents are very robust, in contrast to the cases of two 
components. This happens because the number of com- 
ponents, three, coincides with the number of vertices of 
a triangle. The implication to a colorful vortex lattice 
in the color superconductor is suggestive. Since both the 
number of color and the flavors (up, down, strange) are 
three, it implies a very robust colorful vortex lattice. 

Because of our motivation to simulate a color supercon- 
ductor, we restrict to the case of the symmetric couplings 
in Eq. (pi). Apart from this motivation, we could change 
couplings asymmetrically. If we change either g^ , gij , rrii 
or fii, circulations of vortices are not 1/3 quantized any- 
more. When we consider multi-component BECs as a 
mixture of different atoms, this happens in general. The 
vortex phase diagram in such a case has been studied 
recently for two-component BECs [TS]. These remain as 
future problem. 

If we set 323 = 531 = 0, the third component ^3 decou- 
ples from the others, reducing to two-component BECs 
of ^Pi and 'i'2- There, we have the phase diagram of the 
two component BECs |13ffT5] . It is an interesting future 
problem how the vortex phase of two-component BECs 



is connected to that of a three-component BEC found in 
this paper by gradually increasing the coupling to ^3. 

Another interesting problem is to introduce inter- 
nal coherent couplings between different components via 
Rabi oscillations. In the case of two-component BECs, 
two different fractional vortices are combined by a sine- 
Gordon domain wall to become a two- vortex molecule or 
a dimer [20 . The effects of this term in the vortex phase 
diagram has been studied recently by the present authors 
[19 . For three-components, the introduction of internal 
coherent couplings results in a three-vortex molecule, a 
trimer [23J, and the same for iV-component results in a 
vortex A^-omer [24,. 

Before closing the paper, let us make a comment on 
a possible realization in experiments. Two-component 
BECs of different hyperfine states of the same atom have 
been already realized using the |1,— l) and |2, l) states 
gg and the [2, l) and [2, 2) states gS] of ^^Rb, respec- 
tively. Systems with three components can be realized 
in principle by using a mixture of the above mentioned 
states of ®^Rb via an optical trap [44 , and our prediction 
is testable in laboratory experiments. 
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